Imperial/TP/95-96/3 
|hep-ph/95 102981 
15th October, 1995 

THE CONDENSED MATTER LIMIT OF 
RELATIVISTIC QFT * 



T.S. EVANS"^ 
Theoretical Physics, Blackett Lab., Imperial College, 
Prince Consort Road, London, SW7 2BZ, U.K. 



ABSTRACT 

I study how to apply relativistic quantum field theory to condensed matter systems. 
The motivation for this is examined and then two separate elements are considered. 
First we identify the precise relationship between relativistic and non-relativistic 
fields. Second we consider the need for a chemical potential and how one includes 
this in static and dynamical calculations. 



1. Introduction 

There has been considerable recent interest in condensed matter systems amongst cosmologistsE 
This is because it seems possible to apply methods used in the early universe to some con- 
densed matter systems. The latter calculations can be compared to laboratory experiments 
which in turn can give us further confidence in the early universe results. The table below 
illustrates some of the links which I will discuss in the paper. In particular we note the 
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Table 1. Analogies between Relativistic and Condensed Matter systems. 

theoretical work on cosmic strings, and the links made0 between this and the experimentsi 
on vortices on liquid helium, both ^He and ^He. For example there is the recent work using 
out-of-equilibrium quantum field theory on the creation of cosmic stringg3. In this a single 
complex scalar field was used and it was important to identify the correct physical dynamical 
behaviors, especially the mass, to be used in the calculations. It would be interesting to look 
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at the condensed matter limit as this can be considered as a toy model for liquid ^He. However 
we need to identify precisely what are the non-relativistic fields in terms of the relativistic 
fields. 

First we must specify what we meant by the relativistic and condensed matter limits. 
The energy scales involved in the relativistic problems are such that a particle with typical 
thermal kinetic energy is moving at relativistic speeds, so T>m. Then the Hamiltonian used 
is built around free particles with the usual dispersion relation uj{k) = \/{}? + m^). In more 
practical terms, this usually means that thermal fiuctuations are of order the energy scale 
of the theory, so this is where there is likely to be interesting physics e.g. phase transitions. 
In many cosmological problems, the density of relativistic particles, and relativistic anti- 
particles, n_, is large but essentially the same. In equilibrium we have 

^^^^^ = exp{/3(4):F/.)}' 
so this means that we have a tiny chemical potential. 



n+ — n_ 



<1 =^ M^l (2) 



m 



For instance at high temperatures in the early universe the ratio of the baryon/anti-baryon 
asymmetry to the baryon density is around 10~^°. 

The condensed matter limit is a non-relativistic limit so that T <^ m and for free particles 
the dispersion relation is a;(A;) ~ m + k"^/ (2m). For instance in liquid Helium we are working 
at T ~ IK or less. However, the relevant particles are the Helium atoms whose masses are 
equivalent to lO^^K. The other aspect of the condensed matter limit, and a further difference 
from typical relativistic problems, is that there are essentially no anti-particles so that 

^«l ^ ,^rn. (3) 

It is for this reason that condensed matter texts use a different definition for chemical poten- 
tials. Here I will denote this alternative by Uni- The relationship to the chemical potential, 
II, normally used in relativistic texts, is given by 

Atnr - m » —m. (4) 

Note that for Fermi systems, ii^r has a direct physical meaning in terms of the Fermi 
surface as it is the Fermi energy at zero temperature and plays a similar role at non-zero 
temperatures. For Bosons with no condensation present, is strictly less than the physical 
mass m and it has no direct physical interpretation in terms of some 'Bose-surface'. With 
Bose-Einstein condensation, = m exactly and /i is no longer even a free parameter. The 
treatment of this case is more intricate. I will focus on Bosonic systems in this talk. 



2. Extracting the non-relativistic fields 



In looking at relativistic and non-relativistic bosonic systems, we note that the usual 
equations of motion for the free fields differ by being second order and first order in time 
derivatives respectively; 

(92 + fc2 + m2)$ = 0, (5) 
H5* + ^-M* = 0. (6) 

Note (^) is the Gross-Pitaevskii equation rather than the time-dependent Ginzburg-Landau 
equationi. The latter has no factor of i in front of the time derivative but is generally thought 
by condensed matter theorists to be a better model for systems like ^He. This suggests that 
in trying to express the condensed matter fields \E' in terms of the relativistic fields we must 
split up the relativistic Klein-Gordon equation into two first order (in time) equations. This is 
sometimes done when solving differential equations numerically, when one would take $ and 
9t$ as independent variables. So we will work with $ and its conjugate variable 11 ('~ 
and start not from the usual Lagrangian path integral but from the Hamiltonian version. 
In this way, 11 remains available as an independent variable which we can use, whereas it is 
integrated out in reaching the Lagrangian formulation. Thus a suitable generating functional 
to consider is 

Z[j,f-fi] = J^VU*VUV<^>*V<i>exp(zSe[U,<!>;fi] + zJ{f<l> + j<!>*)y (7) 

where 

B ^ {m = Ht-m (8) 

5e[n, $; /x] = ^ dt J £x (Udt^ + ITdt^*) ~{H~ /iQ) , (9) 
H = J d^x[U*Il+{\/^*){V^)+m^^*^ + 0{\)] (10) 
Q = i J d^x {^*U* -U^). (11) 

and 0(A) represents all cubic and higher order terms. It is useful to think of Se as a matrix 

Clearly proceeding straight from here produces the usual relativistic formalism. So what 
we have to do is to replace the four fields, $, 11 and their conjugates, by four new ones. We do 
this by finding a canonical transformation which leaves us with a temporal first order action 
and which is diagonal in the new fields up to and including the quadratic terms. This is 
sufficient to specify the new combinations to be 

* = \/?^ + ^n* (13) 

V ^ \/2u! 



^ = /|r + -i=n (14) 



together with their complex conjugates. This linear transformation also means that the 
measure in the path integral does not pick up any strange terms 



VIl*VUV<!>*V<^ = D^^P^P^*!)^. (15) 
Physically what has been achieved is easy to see in terms of free fields as we have that 

^free ~ /ci^fca(A?)e~'^-"+'"*, f^free ~ / rf^A? &(^)e~'^"^'"* (16) 

where a{k),b{k) annihilate particles and anti-particles respectively. 

Putting this all together we find that the generating functional now looks like 

Z[J,3*■,^I] = y'p^*P^P^Wexp{S^,„,,e} (17) 

S^,,nr,e = So.nr , /i, ^] + 5o,nr ['J'* , -/i, C^] 

+ (18) 

^0,nr[^*,^,/i,^] = j^'^ dt j (fx^* 



■9 



^ (19) 



on removing total derivatives and in the last expression doing a small momentum expansion. 
Thus the quadratic part of the non-relativistic fields \E' and 'I' has a global 0(2) ® 0(2) sym- 
metry. Particles and anti-particles are separately conserved under this symmetry. However 
typical interaction terms, contained in Sj^^^-, break this symmetry, e.g. 

A|$*$|2 = A^($*^ + vl/*vl/ + vfr*^ + ^*i|r). (20) 

Note that the global 0(2) (g) 0(2) symmetry of non-interacting \E', ^ fields is a mixture of 
the original Poincare symmetries and the original global 0(2) of the $ fields. Put another 
way, the 0(2) ® 0(2) symmetry transformations written in terms of the $ fields are non-local 
transformations. If we make the transformation 

^ ^ 6^^^+")^, ^ e^^''-^)^ (21) 

then this is equivalent to making the transformation 





(22) 



An alternative way of looking at this approximate 0(2) C?) 0(2) symmetry is to make a 
transformation to \E'', ^' defined in the same way as before except that we replace uj by the 
in (p!3|) and (p^). This version is mentioned by Haber and WeldonQ (see the sixth 



mass m 



citation/footnote). Then these ^' fields show only an approximate 0(2) (8)0(2) symmetry, 
the symmetry being broken by terms of order k/m as well as by the interactions in Si but 
the symmetry is now a standard global symmetry in terms of the $ and 11 fields. Since 
the transformation from $ to \E'' is global, this approximate 0(2) ® 0(2) symmetry is also 
a approximate global symmetry of the original $ Lagrangian. It is not usually discussed 
because the terms which break the symmetry are large for relativistic particles. 

So far everything has been exact, for any value of fi and T. In practice for a condensed 
matter system m c:^ fi ^ T. The low temperature means that we can make the non-relativistic 
approximation 

r-iP r [ d 1 

5o,nr[^*,^,/i,^] ^ dt J d^x<il*^i— + —V^ + ^-m 



^ (23) 



The high particle density, or equivalently the large chemical potential, means that the mass 
gap term in the anti-particle ^ free Hamiltonian dominates as m — fir^T'^rrir^fi. This 
means that the quadratic part of the anti-particle Hamiltonian is to a good approximation 
2m^*^. On integrating out the anti-particle fields we find that we just set = ^ = and 
so the anti-particle fields just drop out. 
This then leaves us with 



-il3 



■9 1 „2 



^ + (24) 



where for A|$*$p we have S'/[\E'*, \['] = A(\E'*\l')^/(4m^). This means that we are now left with 
an approximate theory, (p^, which also has an exact 0(2) symmetry. In terms of the original 
theory, this new 0(2) is broken by the terms of 0{k/m),0{{fi — m)/m) which are dropped in 
reaching (|2^). It is a different symmetry from the original exact global 0(2) of the relativistic 
theory. The latter mixed particles and anti-particles, whereas (^3|) does not. 

There is one caveat. The <I>-fields are periodic in imaginary time, but the Il-fields are not. 
This means that the ^'-fields can only be approximately periodic, although when (^ip is used 
exact periodicity is assumed. 

The case of fermionic fields is rather different from this bosonic case as the Dirac equation 
is already first order in time. The problem with fermions is rather that one must generate the 
V^/2m term from the first order 7.V term, and also go from four by four Dirac gamma algebra 
to the two by two Pauli matrices. Essentially in this case one must follow a Foldy-Wouthuysen 
approach, but this will be considered elsewhere. 

3. Chemical Potential 

There are two approaches to the chemical potential in thermal field theory, which are 
distinguished by the way they split the physical information between initial conditions and 
microscopic dynamics - see table H. 





Microscopic Dynamics 


Initial Conditions 


Method I 






Method II 


H 


Hto) = e^^'<I>(to - tP) 



Table 2. Two Approaches to Chemical Potentials. 

In Method I, the H and /iQ terms of the density matrix are merged into a single an 
effective Hamiltonian, He which then explicitly depends on the chemical potential 

p = e-^(^-'^«) = e-^^^^ (25) 

One of the great advantages of this method is that this new effective theory looks just like a 
theory with no chemical potential. This makes the boundary conditions the usual simple (anti- 
)periodic ones for bosonic (fermionic) fields so that all the usual well established methods can 
be called upon. Not surprisingly this approach is very popularS'ii. While the initial picture is 
simpler, it is often the case that new complications due to the presence of a chemical potential 
arise further down the road. All fields, and hence all Green functions, in Method I have a 
subscript e to denote that is being used. 

In Method II, all aspects of the equilibrium density matrix are treated in the same way. 
Thus one encodes both the temperature and the chemical potential in the boundary con- 
ditions, and the Hamiltonian used is physical one with no information about the density 



matrix included. This leads to a rather more complicated boundary condition, which for 
bosonic fields is of the form 

<l>(to) = e^'^$(to-^/5). (26) 

This makes matters more difficult than the zero chemical potential case from the very start. 
The advantage of this approach is that the density matrix is often best thought of as defining 
the initial conditions. The microscopic dynamics of the physical fields, as described by H, 
has nothing to do with the initial conditions. Thus keeping the Hamiltonian used free of 
contamination from the initial conditions means that Method II is £Lvery physical approach. 
This is the approach used in the review of Landsman and van Weert". It is also interesting to 
note that it effectively the way that standard condensed matter textslllll deal with chemical 
potentials associated with fermionic degrees of freedom. 

In fact for most out-of-equilibrium problems, the reverse is also true, that is the initial 
conditions need have nothing to do with the microscopic dynamics. Indeed in Closed Time 
Path approaches to non-equilibrium problemsHlii^ the initial conditions may always be en- 
coded in terms of some 'temperature' and some 'Hamiltonian' as we can always write the 
density matrix as p = exp{— /3j„jtifj„jt} . However there is no reason why this Hinu has any- 
thing to do the Hamiltonian which controls the microscopic interactions and hence evolution 
of the system. 

Of course in many practical situations the initial conditions are related to some physical 
Hamiltonian. Equilibrium is an extreme case where the same H appears in initial density 
matrix and evolution. However, Method I is still using a mathematical trick to 'simplify' 
matters but it is doing it at the expense of physical clarity. The discussion above suggests 
that problems with Method I may occur when looking at dynamical problems. This is because 
H not Hf. really controls the microscopic evolution yet in Method I Hf. is used for everything. 
We will now justify this assertion in detail by looking first at the propagators which appear 
in both methods, i.e. we will solve simplest Droblem, that of a single free complex scalar field. 

We will start by following Method I@'EE2l, which is to use real fields 0ei,0e2, where <l>e = 
(0ei + ^0e2)/y2. Thus the action, S'e, is of the form 

(27) 




'^'-/i^el) -^(V0e2)'-^mVe'2 + O(A) 



2 \ dt 2' ^ ' 2 



dtdt' {(l)el{t)Ae2{t)),Dab{t;t') ^ +0(A). (28) 

The inverse matrix propagator, D^^, is given in energy-momentum coordinates, less an overall 
conservation delta function, by 

^ 2..„/,3. k) ^ ( ) (29) 

where for notational simplicity we have used the imaginary time form for the propagator. This 
propagator is not diagonal in the one and two fields but it is easy to find the two eigenvalues 
of D^^ which are ffco + /i)^ — oj"^ and (fco — /i)^ — w^. This suggests that there are four poles 
in the propagatoru at /cq = ±L(J ± /z. However, we started with a single complex scalar field so 
we expect only two poles. What has happened? 



The solution is to look more carefully at the diagonalised Lagrangian. What we find is 
that the eigenstates are (1,«) and and hence the diagonalisation process has merely 



moved us from real fields 



')e2 back to the original complex fields $e,$*. This is not 



surprising as we know that the latter are eigenstates of the charge operator, the real fields 
are not. It makes more sense in this case (where the symmetry is unbroken) to stick with the 
Q eigenstates. If we write out the free action after we have diagonalised (ESl), we obtain 
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U,me\t,t')<i>:{t' 
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dtdt' <^l{t)G-\t,t')<^e{t') 



PV JO 

where the two propagators G and H are defined and related through 

tGe{t,t') = Tr{e^(^-^^T[<I>,(t)<I):(t')]} 
tH,it,t') = Tr{e^(^-'^'?T[$:(t)$e(t')]} = ^G'e(t',t) 
G;\ko,k) = [{ko + fif-u;'] = H;\-ko,k) 



(30) 
(31) 
(32) 



(33) 
(34) 
(35) 



were all the times are time ordering are Euclidean valued, and the energies are taken at 
ko = 27m /f3 where n is an integer. Thus we see from ( |30|) that what appears to be four poles 
is actually the same two poles repeated a second time but with the time or energy argument 

521) much quicker by working with $, $* from 



reversed, as (^) shows. In fact we can obtain 
the start. This should not be surprising as these are the fields which are the eigenstates of 
the charge operator. Anyway there are only two distinct poles as one expects given that we 
have a single complex scalar field. For Ge, (pSI), the poles are located at 



kn 



(36) 



Now let us turn to Method II and compare the propagator found there with (|55|). Essen- 
tially all we have to do is solve the Klein-Gordon equation subject to the weighted periodic 
boundary condition (^Bp. Of course one needs other boundary conditions but these are pro- 
vided by the equal time commutation relations. It is simple to look up the answerO and the 
time-ordered two-point Green function, Gt, is given by 



iGt(t, X] t' , x' 
1 



(2vr) 



Trie^'^^-^^rmt, x)^*{t', x')]} 



d'^k eyi'i){—ikot + ik.x} 



kn 



2TT6{ko — uj) 



2n6{ko + uj) 



exp{l3{ko — /i)} — 1 2uj 



exp{f3{-ko - /i)} 



2u 



(37) 



Thus we see clearly that in Method II there are only two poles and they are at the same place 
as they are for a free field at zero density and arbitrary temperature, namely at 



ko = ±uj. 



{31 



So in Method II we find that the poles of the free propagator do not depend on any of the 
many-body parameters, and are only dependent on the form of the quadratic part of the 



Hamiltonian. Thus the position of the poles depends only on microscopic physics which is 
described by the true physical Hamiltonian. On the other hand, the weighting or residue of 
the poles does depend on the initial conditions. Here the real-time propagator shows us this 
physics very clearly. We have a zero temperature term which is a pure quantum fluctuation 
contribution so has no explicit /i or /3. Then there is one particle and one anti-particle 
pole, each weighted by the appropriate statistical distribution. These represent statistical 
fluctuations due to the real particles in the heat bath. This is why the particle pole has a 
much larger residue than the anti-particle pole when there are many particles but few anti- 
particles (large positive /x). So Method H leaves us with a picture which fits our physical 
intuition. 



However, we now have an apparent contradiction between the results of Method I (|3q ) 
and Method H (|38| ). As the investigation of the relationships between various Green functions 
and the Green functions calculated in different formalismst3 showed, it is imperative that we 
are extremely careful in identifying the Green functions we are studying. 

It is immediately obvious that the difference encountered here between Methods I and II 
is not a result of looking at different types Green functions. In the analysis of Method I we 
looked at Euclidean two-point Green functions and, implicitly, their analytic continuations, 
the retarded and advanced propagators. It is well known though that the real parts of these 
two-point functions are identical to the rea]jia|diS of the time-ordered Green function, which 
was considered when looking at Method IiyMyi. So we must go back to basics to find the 
source of the contradiction. 

As we are looking at a propagator in Method I, let us write down the generating functional 

Ze[3l,3e] = Tr{e-^(^-^<3)2-[e-^/{i:*e+**>)]| ^39) 
= 5: ^,,($,to-^/5|re-^^(^^*^+*^>)|$,to) (40) 

where we have used the unphysical effective Hamiltonian Hf. = H — nQ of Method I to evolve 
the bra. This means that in order to calculate Zg, say by using path integrals, we have to 
evolve all fields and states using ifg not the physical H. Thus all the fields which appear in 
Zg, or its derivatives, satisfy in the Heisenberg picture 

$^(t) = e^^=*<l>(0)e-^^=*. (41) 

The fields appearing in Method I are therefore not the same for an arbitrary time as the 
physical fields $, which satisfy in the Heisenberg picture 

$(t) = e^^*$(0)e-^^*. (42) 

This is the reason we have been carefully labeling the fields used in Method I with the e 
subscript. Likewise, the Green functions must be distinguished. 

There is however a simple link between the two fields, and we see that from ( ^if) 

$^(t) = e^(^-^'3)*<l>(0)e-*(^=-''^)* (43) 

= e-^^^*$(t)e^^^* (45) 
^^{t) = e*^'^*<I>(t). (46) 



where we have assumed that our field has a charge of q. The same connection is also noted 
elsewhereEl but there is no discussion of Method II as a practical calculational tool for bosons. 



Now, to make a link between the Green functions, we need to make a connection between 
the generating functionals in Method I and II. So we choose the sources such that 



(47) 
(48) 



Jeit) = e-'^'j{t). 

as this then gives 

z[f,j] = z4j:,j,]. 

The generating functional in Method I, Z^, was given in ( |50| ) and that of Method II is 

Z[j*,j] = Tre-^(^-'^«)Te-'/(^**+**^) (49) 
= H{^,to-tP\Te-^I^'"'-"'*^%to) (50) 

This means that Method II generates physical Green functions e.g. 



r(2^)(ii,---,^iv;si,...,sjv) 

= Tr 



dj*{ti) . . . dj*{tn)dj{si) . . . dj{sN) 



(51) 
(52) 



where the T indicates path ordering appropriate for whatever finite temperature formalism 
is being considered. Likewise we see that with Method I we generate thermal expectation 
values but of the unphysical fields ^e,^*e, namely 



di2N)z 



dj*{ti) . . . dj*{tn)dj{si) . . . dje{sN) 



Tr 



(53) 
(54) 



Now we can use (^), (^3) ^'^'^ (HD show how the unphysical Green functions of Method 
I, Fe of (0), are linked to the physical ones of Method II, V of (^21). We find 



■^«)T[<|.e(tl)...$e(t7v)$:(Sl)...$e(Siv)] 



r(2^)(ti,...,t^;.i,...,.;v) = e-^^S.(*^-^^)r('^)(ti,...,t;v;«i, 

In particular, we find that for the propagators we have that 



Sat) 



D{t,s) 



D{h) 



Tr 



A'Q)T[$(t)<l>*(s)] 



De{ko - /i) 



(55) 

(56) 

(57) 
(58) 



These relations apply to any approximation to the full Green functions and in particular for 
the free field case. This resolves the paradox we found above, (^6[) vs. (|3|). Equation (|5^ ) 
tells us that the poles in Method I and Method II are related by a simple shift of fi in the 
energy. 

It is important to note from the relation between the two generating functions that the 
partition function is identical in both cases 



Z = Z[0,0] = Ze[0,0] 



(59) 



Thus both methods give the same answer for the partition function and so will give answers 
directly for all thermodynamic quantities. Indeed, any static quantity, i.e. time independent 



or equivalently zero energy quantity, can be extracted directly from either method. Only 
when one has time dependence or is looking at non-zero energy does one need to use formulae 
such as (^5]) to move between the two formalisms. 

Finally we mention a complication with symnietrv breaking, or equivalently Bose-Einstein 
condensation. This area has long been of interesta'Qfl'clll^ and has had a lot of recent attention 
ult of the experimental observation of Bose-Einstein condensation in "Quantum Atom" 
systemsEj. So far we have been working in a situation where the symmetry is unbroken, 
($) = 0. However what happens if we are in a phase where this is not true? In Method I 
with periodic boundary conditions we proceed as usual and shift the fields 

$e(t,x) = Ve + ^eit,x) (60) 
= V, (^e) = 0, (61) 

where Ve is a constant in the usual way. However in Method H such a simple approach is not 



compatible with the more complicated boundary conditions (|2^). By using the link between 
the fields in Methods I and II, (|46|) we see that we have to use 



$(t,f) = f(t) + *(t,f) (62) 
($e)=t;(t), vit) = e'^^'ve, (^e) = 0. (63) 



4. Conclusions 

We have seen how to link the fields used in relativistic and non-relativistic theories. Equa- 
tion (|4^) gives this link. Note that we worked within Method I but the formula works for 
Method II as weU. 

I have also discussed how to link the two approaches used when including chemical poten- 
tials. In Method I, fields <I>e are used which evolve according to = H — ^Q. This means 
that for these fields, we are measuring the energy on a scale shifted by /i relative to the stan- 
dard relativistic field theory definition. Thus Method I is actually quite close to Method II. 
Method I is only 'unphysical' in the sense that it is measuring its energies from a zero which 
is shifted from the standard of relativistic field theory,^ Method I is favoured by condensed 
matter physicists when dealing with lioaDnic systemgl3. It is also used in relativistic texts 
for calculations of partition functions@'clliy. It is easy to use for free energies because of the 
equality (^Up between the partition function of Method I, and that of method II. 

However, for dynamical quantities there is a difference between Methods I and II, sum- 
marised by (|55|). In this case. Method II is clearer in my opinion as it sticks rigidly to the 
conventions of relativistic field theory. In particular all energies are measured relative to the 
same zero, which is the energy of the vacuum state (in a simple unbroken case). However 
the same physical information can be obtained from Method II. The problem is that energies 
are then measured relative to a zero shifted by from that of Method I and the standard 
relativistic zero. This means that fields with different charges under the various symmetry 
groups involved will have their zero's of energy shifted by different amounts! Thus in try- 
ing to understand the electroweak model one runs a serious risk of getting confused about 
the energies involved of various modes. So in discussions of dynamics, poles of propagators 
or whatever, in a relativistic field particle context but where Method I is used, such as in 
Bernstein and Dodelsona, one must be very careful that one understands that the energy zero 
has been shifted. Shifting your zero of energy does not mean that anything interesting has 
happened physically! 



However it is a matter of personal taste which method one chooses to use. The same 
physical information is contained in both approaches. One must just be very careful when 
interpreting the information and be very clear what energy zero one is working with. Method 
I is closest to the way condensed matter physicists work with bosonic systems, with energies 
measured relative to the Fermi surface (which is equal to // at zero temperature) . Method II 
means that we stick to the conventions used amongst relativistic particle physicists. 
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